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Strong and weak gravitational field in R + /i 4 / R gravity 
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Abstract We introduce a new approach for investi- 
gating the weak field limit of vacuum field equations 
in f(R) gravity and we find the weak field limit of 
f{R) = R + Li 4 /R gravity. Furthermore, we study the 
strong gravity regime in R + /i 4 / R model of f(R) grav- 
ity. We show the existence of strong gravitational field 
in vacuum for such model. We find out in the limit 
l-i — > , the weak field limit and the strong gravitational 
field can be regarded as a perturbed Schwarzschild met- 
ric. 

Keywords Spherically symmetric solution. f(R) grav- 
ity. General relativity 



1 Introductions 



Observations on supern ova type la (jRiess et al. 111998 



part of Hilbcrt-Einstein lagrangian ( Capozziello 


2002 


Carroll et al. 2004, 


20051 


Clifton and Barrow 1 


2005 


Noiiri and Odintsov 


20031 


Sawicki and Hu 12007b 


Evans et all 20081: As 


timohammadi et al. 2009). 


One 



of the initiative f(R) models supposed to explain the 
positive acceleration of exp anding universe has f(R) 
action as f(R) = R - /u 4 /R dCarroll et al. 112004 . Af- 
ter proposing the f(R) = R — fi A /R model, it was 
appeared this model suffer several problems. In the 
metric formalism, initially Dolgov and Kawasaki dis- 
covered the violent instabil ity in the matter sector 



( Dolgov and Kawasaki 1 120031 ). The analysis o f this in 



IPerlmutter et al. Ill999l ). cosmic microwave bac kground 
(ISpergel et al. 1120031 ) and large scale structure ( Tegmark 



stabi l ity generalized to arbitra ry f(R) models (jFaraoni 
20061: ISawicki and Hu ll2007al) and it was shown than 
an f(R) model is stable if d 2 f/dR 2 > and unsta- 
ble if d 2 f/dR 2 < 0. Thus we can deduce R - ^/R 
suffer the Dolgov-Kawasaki instability but this insta- 
bility removes in the R + [i 4 / R model, where [i 4 > 0. 
Furthermore, one can see in the R ~ [i 4 j 'R model the 



e -(- cgsraology is inconsistent with observation when non- 



2004]), all indicate that the expansion of the universe 
is not proceeding as predicted by general relativity, if 
the universe is homogeneous, spatially flat, and filled 
with relativistic matter. An interesting approach to ex- 
plain the positive acceleration of the universe is f(R) 
theories of gravity which generalize the geometrical 
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relauvistic mat ter is present. In fact the re i s no matter 
domi nant era ( Amendola et al. 2007al lb: Evans et~al 



20081) . However, the recent study shows the standard 
epoch of matter domination can be obtained in the 
R + H 4 /R model ((Evans et alll2008l) . 

It is obvious that a viable theory of gravity must 
have the correct newtonian limit. Indeed a viable the- 
ory of f(R) gravity must pass solar system tests. After 
the R — [i A /R was suggested as the solution of cosmic- 
acceleration puzzle, it has been argued t hat this theor y 
is inconsistent with solar system tests (jChiba 1 l2003l) . 
This claim was based on the fact that metric f(R) grav- 
ity is equivalent to u> — Brans-Dickc theory, while 
the observational constraint is uj > 40000. But this 
is not quite the case and it is possible to investigate 
the spherical symmetric solutions of f(R) gravity with- 
out invoking t he equivalence of f(R) grav i ty and scalar 
tensor theory ( Clifton and Barrow I l2005i ICembranos 



2006 


Sawicki and Hu 


2007ti Multamaki and Vilia 


The analysis of spherically symmetric solution can be 
carried out using schwarzschild coordinate 


2006 


CaDOzziello et al. 


20081 120091: Saaidi et al 12010: 


Capozziello et al. 2010 


). It has been shown that 



space-ti me as a spherical symmetric solutions of field 
cquation (|Multamaki and Vilia 1 120061) . Hence R— / 1 R 



model has a Schwarzschild-de Sitter solution with con- 
stant curvature as R = where this is not the case 
111 R + n A /R model. 

In this paper we study the R + fJ. 4 /R model of f(R) 
gravity. We find the static spherically symmetric solu- 
tion of vacuum field equation in both weak field limit 
and strong gravity regime, moreover, the weak field 
analysis can be expanded on f(R) models of the form 
f(R) =R + eh(R). 



2 Weak field limit 

In this section we investigate the weak field solution of 
vacuum field equation in f(R) theories of gravity. We 
are interested in model of the form f(R) = R + eh(R), 
with e an adjustable small parameter. The motivation 
for discussing these models is that the nonlinear cur- 
vature terms that grow at low curvature can lead to 
the late time positive acceleration, but during the stan- 
dard matter dominated epoch, where the curvature is 
assumed to be relatively high, could have a negligible 
effect. 

The field equations for these models are 



G 



R m\ 

<p(R)J 



<p(R) + kT^, 



where ip(R) = dh(R)/dR. Contracting the field equa- 
tion we obtain 



R 



R 



2h(R) 
<p(R) 



3D 



ip{R) - kT. 



(2) 



Where for the vacuum T^ V ,T = 0. If e = the above 
equations reduce to Einstein equation. Hence we sup- 
pose G^y and R in the r.h.s of Eqs. (|ll2[) can be ne- 
glected for small values of e. Furthermore if the con- 
dition lim [h(R)/<p(R)] — is satisfied we can neglect 

this term too. Neglecting these terms leads to the fol- 
lowing equations 

G>„ = -c [g^U - V^V,] <p(R), (3) 
and 

R = e3Uip(R). (4) 



In the weak field limit approximation the metric de- 
viates slightly from the Minkowski metric, so we can 
write 



A(r) = l + a(r), 
B(r) = l + 6(r), 
|a|,|6|«l- 



(6) 



When solving the field equations (|3l4p we will keep only 
terms linear in the perturbations a(r), b(r). Hence 
equations (|3l4p leads to 



a' 

— + 


b 


r 


^2 


V 


b 


- + 

r 




and 





_ 2d<p{R) 
r dr 



(7) 



R = 3eW 2 (f(R). (8) 
where (/) indicates a derivation with respect to r. 
2.1 f(R) = R 1+£ 



This model is considered in ( Clifton and Barrow 120051) . 
It is shown that this model has an exact spherically 
symmetric vacuum solution and regarding the general 
line-element in Eq.([5]), it may be written as 

(1) M r \ = r 2e(l+2e)/(l-e) +c r -(l-4e)/(l-e) ; 



B(r) 



(1 



(1 - 2e + 4e 2 )(l - 2e - 2e 2 ) 
x ^l + cr -(i-2 £ +4 £ 2 )/(i- e )^ 



where c is a constant. In the limit e — > 0, these solutions 
become 



ds 2 



l + 2elnr+- ) dt z + 1 + 2e + - ] dr 
r J V r 

2 ir\2 



^dVL\ 



(9) 



because we seek the weak field limit, in above equation 
we assume c/r <C 1. 

Since we are interested in the limit e — > 0, we may 
expand f(R) = R 1+£ around e = 0. Then we have 

f(R) = R + eRlnR, 
h(R) = R]nR, 
<p(R) = 1 + lni?. 
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It is clear that h(R) satisfies the condition 
lim \h(R)/<p(R)] = 0. 

ii->0 

Inserting (p(R) in the trace equation ([5]), the Ricci scalar 
is obtained as 



R 



6e 



Then we arrive at the solutions of Eq. ([7]) 



a = — h 2e In r, 

r 



b=-+2e, 
r 



(10) 



(11) 



where c is a constant. We can see our solutions are in 
agreement with the exact solutions ©■ Also one can 
check neglecting R, G M „ and h(R)/tp(R) in Eq.([IJ[2]) is 
reasonable. 

2.2 f(R) =i? + elni? 

For this model f(R) = 1/R. Solving trace equation 
and field equations ([7]) we obtain 



Indeed by this approach the Post-Newtonian parame- 
ter is found as jppn = 1/2 while the measurements 
indic ate ^ PPN = I + (2.1 ±2.3) x 10~ 5 jBertotti et al. 
2003). Also we must note that using equivalence be 



tween f(R) gravity and scalar tensor gravity one can 
find models which are consistent with the solar sys- 
tem tests. This consistency can be made by giving the 
scalar a high mass or exploiting the so-called chameleon 



effect (IMota and Barrow 1 12004c iKhoury and Weltman 
2004; ICapozziello and Tsujikawa l2008t iFaulkner et al. 



20071) . However, when one is using equivalence between 
f{R) gravity and scalar tensor gravity, the continuity 
of scalar field or its equivalent, the Ricci scalar, at the 
matter boundary is crucial condition which is not the 
case in Einstein gravity. But in this work we don't 
adopt the continuity of Ricci scalar for solving the field 
equations. Instead, we suppose that when /i tends to 
zero we arrive at the Einstein gravity. Thus we find a 
solution for 1/R m odel which is radically different from 
other solutions in ( Erickcek et al. 2006t Chiba et al, 
20071) . 

For this model we have 



R 



and 



6e 



a = b = — 



2M 



(12) 



(13) 



where M is a constant. Therefore the space time metric 
for empty space in this model is 



ds 2 



1 



1 



2M 
r 

2M 
r 



-r | 



dr 



(14) 



We can see, the generalized Newtonian potential is 
M 1 



4> 



G 



(15) 



This generalized gravitational potential has two terms. 
The first term is the standard Newtonian potential 
and the second term make a constant acceleration, 
+ ^/e/24, which is indepen dent of the mass of star. In 
( Saffari and Rahvar 1 120081 ) this metric is used to ad- 
dress the Pioneer's anomalous. 

2.3 f(R) = R± f i i /R 

Based on equivalence between f(R) gravity and Brans- 
Dicke theory with uj = 0, it was argued t hat this theor y 
is inconsistent with solar system tests (jChiba I I2003T ). 



h(R) = ±1/R, 
ip(R) = Tl/R 2 , 



(16) 



where h{R) fulfills the condition Urn {h(R)/ip(R)} = 0. 
Solving Eqs. (|7l8[) we obtain 



R = =F7a/i 3 r 3 ; 

U 4 1 4 4 

i? 2 49a 2 ^ ' 

21 , 3 4 4 

a = ± —afi^r 3 , 

r 4 

2M 4 4 

o = ± u[i 3 r 3 . 

r 



(17) 



where a 3 = 4/147 and M is a constant. Therefore the 
metric for space time is 



ds 2 



, 2M 3 4 4. 

1 ± -«ii 3 r 3 dt 

r 4 



2M 



1- — ±a^ir« ) dr 2 +r z dVL\ (18) 
r 



We can use the isotropic form, by introducing 
a new radius, p, which defined as 



I 2M , 3 4 4 
r = pJl + — ± -a^3 p3 , 

and therefore the equivalent metric can be read 
as 

, 9 / 2M 3 4 4 , 9 

ds 2 = -(1 - — — ± -ap3 p 3)dt 2 



■1 



+ 1 + ±-a^pt )(d P 2 + P 2 dn 2 ). 

P 4 

From above metric one can see in the asymptotic 
behavior, p — > 0, 7ppw ~ 1 can be obtained. 
From Eq. (fTT|) it is obvious that in the limit p — > 0, 
p A /R 2 tends to zero, so there is not singularity in the 
field equations. Also one can check neglecting R, G^, 
and h{R)/ip(R) in Eq.([TJ[2]) is reasonable. 

2.4 Interior solution in the f(R) = R + p 4 /R model 

In this section we discuss the interior gravitational field 
in the spherically symmetric case of static mass distri- 
bution in the f(R) = R + p A / R model where p — > 0. 
So we seek a spherically symmetric, static solution and 
we adopt the metric©. In this model we may rewrite 
field equation ([1]) and trace equation @ as 

G; = (S'R + G' + S-n-V^^ + kTZ, (19) 



where in the above equation we used the empty space 
solution Eq. (fl"8]) . From the above equation we may 
determine the parameter M. It is seen that in the 
p — > limit this constant reduces to the Schwarzschild 
radius. Furtherm ore, according to cos mological studies 
p 2 = 10~ 52 m -2 (jCarroll et al. 1 120041 ) so, regarding a 



typical solar system, in Eq. (l26[) we may neglect terms 
at order O (ppfr j ■ 

From equation (|26|) it is obvious that the physical 
interpretation of the parameter M differ from that of 
general relativity. Also from the above equation it is 
clear that in the 1/R gravity the external solution de- 
pends on the shape of matter distribution. 



3 Strong Gravity Regime in R + £t 4 /-R Model 

In this section we investigate the existence of strong 
gravitational field for f(R) = R + p 4 /R model of f(R) 
gravity. We can rewrite the field equation ([T]) as 



R = 3(R-a)^--kT. 

R z 



(20) 



From Eq.© it is obvious that as p — > 0, R — > — kT, so 
assuming p A <C —kT, in the r.h.s of Eq. (|19[) we may ne- 
glect those terms that contain p 4 /R 2 . Thus field equa- 
tions (fT9"f reduce to Einstein equations hence we may 
write 



(21) 



furthermore the conservation equation, T^ v = 0, 
leads to 



A! 



(22) 



where p, p are pressure and density of matter. To obtain 
metric components (A, B), we use Eq.([22|) and rr and 
tt components of Eq. ([2~1"T) 



A' B B-l , 
Gl = - 1 — ~ hp, 



Gr 



A r r 2 
B' B-l 



-kpc 



r r 
Solving Eq.(|23p we may write 

1 f r 
B = 1 - -kc 2 / p(x)x 2 dx + O 
r Jo 



( P 4 



V fc 2 T 2 



(23) 
(24) 

(25) 



From continuity of the metric component B(r), on 
the boundary surface r = r$ we find 

— jf p{x)x 2 dx + ap3r i + O f^j - — , (26) 



^ ( R~ 2 



where we have used the trace equation 
R= -3[R+a] (p 4 /R 2 ) . 



(27) 



(28) 



In the above equation we have neglected the energy- 
momentum tensor of matter because we investigate the 
strong gravitational field around a spherically symmet- 
ric distribution of matter. Adopting the general spher- 
ically symmetric metric (|5|), we can rewrite the trace 
equation (|28|) and (rr),(tt) components of field equa- 
tion (H3 as 



dr 2 

(pVr 2 ) 



2 d 
r dr 
R 
3' 



BA' 
~A~ 



d 
dr 



R 

(29a) 



,•2 



+ 



1 

B' d 
YJr 



[i-p 4 /R 2 



(BA! M B-1 
\~r~A M 

B & 2 + 



R 



(29b) 



BA' d 
~2A~dH- 



[p 4 /R 2 



R 

"3' 



(29c) 



where(/) denotes derivation with respect to the (r). 
In the previous section we showed, (R + p A / R) model 
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Fig. 1 a against x. The red-solid line shows numerical re- 
sults of Eqs.(32). The green-dashed line represents approx- 
imate solution for x <C 1 (Eq.(34a)) and the black-dotted 
line is the approximate solution for x 3> 1 (Eq.(36a)). A 
close up on the origin of main figure is presented . 



Fig. 2 b against x. The red-solid line shows numerical re- 
sults of Eqs.(32). The green-dashed line represents approx- 
imate solution for x <C 1 (Eq.(34b)) and the black-dotted 
line is the approximate solution for x 3> 1 (Eq.(36b)). A 
close up on the origin of main figure is presented . 



has the week field solution as 
2M 3 



ds 



1 



1 - 



2M 



4 

a(fj,r) 



a(fj,r)~- 



dt 2 



dr 2 +r 2 dn 2 , 



(30) 



where a = (4/147) 1 / 3 . It is obvious this metric re- 
duces to Schwarzschild metric in the limit fx -+ 0. 
Now we seek the solution of field equation in the limit 
(r -+ 2M). Without loss of generality we can assume 
2 M = 1. In order to solve equations (J2U) we use some 
definitions as 

4> = i/R, 



i = -m 4/3 , 

A = l-l 

r 

B = 1 - - +7&(r). 



•7a(r), 



(31) 



Because we seek the solution in the limit r — > 1, we 
may define a new variable as x = r — 1. Using these 
definitions we can rewrite Eq. (|29[) as 



2b 



b' + a' (x + !)(&- a) 



dx 
1 



(.T + l) 2 

( x d 2 



V 1 
- a' 



2 2(x + ja{x + 1)) 

i0 2 _ fl ,o\ 1 
dx \3 

2x + 1 d 



1<P 



1 dx 2 (x + l) 2 dx 



(32a) 



a rf 2 1 d ( d 2 

x + ldx 2 + 2{x + l) 2 dx~ +7 V ffe 2 " 



2"cfa. 



30 + ^(.T + l) 2 

(*TIF + V 



c + 1 a; + 7a(x + 1) 

(1 + 7 2 ) (32b) 



1 / (» + l)(6-o) 

2 V x + ja(x + 1) 

1_ ( b_ 

~ 30 + \ (x + l) 2 + x 



1 



{x + iy 



1 



(l + 7 2 ) 



dx 
(32c) 



where (/) denotes derivation with respect to the (x). 
For the limit fx — > 0, in the above equations we suppose 
that we can neglect terms containing 7 . After solv- 
ing equations we check this assumption. By neglecting 
these terms, equations 1321 can be rewritten as 



1 

30 



1 dx 2 



2x+l d 
(x + l) 2 dx 



(33a) 



b a' b — a 1 

{x + 1) 2 + x + 1 + x(x + l) 2 ~~ ~30 
x d 2 Id 



x+ldx 2 2(x + l) 2 dx 



1 d<t> 2 



b' 



2 (x + 1) 2 dx 30 {x + iy 



x + 1 



(33b) 



(33c) 
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Fig. 3 ip versus x. The red-solid line shows numerical re- 
sults of Eq.(32a). The green-dashed line represents approx- 
imate solution for x <^ 1 (Eq.(34c)) and the black-dotted 
line is the approximate solution for x 3> 1 (Eq.(36c)). A 
close up on the origin of main figure is presented . 



conformal transformation and changing coordinate we 
can see the strong field solution ([35]) is 

ds 2 = 
+rW, 

which is valid in the range of (2M /i) 4 < r/2M - 1< 1 
and farther where r 3> 2M, the metric of space time 
can be approximated by the metric (|30|). Furthermore, 
we have solved field equations (j32J) numerically and pre- 
sented the results in figures [TJ[3J and|31 The plots show 
that the numerical results are in agreement with the 
analytical solutions ()34I36|) in their region of validity. 



In the limit i<1, solutions of Eq. (|3"3"]) are 



3 /4 



00 = 8 U ' X 



1/3 



2/3 



b 



1 fA 



1/3 



,2/3 



1/3 

1 4 1 x 



1/3 



(34a) 
(34b) 
(34c) 



Thus we obtain the metric for i < 1 as 
ds 2 = 



1 



r 2 dft 2 . 



(35) 



Furthermore, for i> 1, we can obtain the solutions of 
equations (f33|) as 



4/3 



boo OLX ^ , 

1 



7a 



x 2 '\ 



(36a) 
(36b) 
(36c) 



which arc in agreement with week field limit (|30|) . Now 
we can check the validity of our assumption. Consid- 
ering the solutions (f3"6"|). shows that neglecting terms 
containing 7 in Eqs. (|32[) is valid only for x 3>| 7 s | 
or 1 » /1 4 . Hence the metric (|34|) is solution of field 
equations in the range of /x 4 <C x <C 1. By performing a 



4 Discussion 

We studied spherically symmetric solution of f(R) 
gravity. At first a new approach for investigating the 
weak field limit of vacuum field equations in f(R) grav- 
ity was introduced. Our results for the weak field 
limit of some studied f(R) models are in agreement 
with their known solutions. We solved the field equa- 
tions for f(R) = R + fJ. 4 /R gravity at weak field limit 
and obtained a solution which differs slightly from the 
schwarzschild metric. Our results are against the ar- 
guments that f(R) models are ill defined because of 
the equivalence of f(R) gravity and Brans-Dickc grav- 
ity with u>bd = which leads to -fppN = 1/2. In fact 
our results are in a greement with the recen t article of 



Capozziello et al. ( Capozziello et al. 201fi ). in which 



they have studied Newtonian limit of the f(R) gravity 
by considering that fourth order gravity models are dy- 
namically equivalent to the O'Hanlon lagrangian and 
they have shown fourth order gravity models can not 
be ruled out only on the base of analogy with Brans- 
Dickc gravity with lubd = 0. Moreover, regarding the 
results for the weak field limit, we investigated the 
strong field regime for this model and showed that if 
(r - 2M)/(2A/) 5 > ^i 4 , where r and 2M arc radius 
and Schwarzschild radius in the Schwarzschild coordi- 
nate respectively, the gravitational field is a perturbed 
Schwarzschild metric even in strong gravity regime, fi- 
nally we solved the master equations numerically by 
setting the initial value conditions using the analytical 
answers of the strong gravity region. In figures (JTJ and 
we plotted the analytical and numerical solutions 
of the components of the metric, a and b, versus ra- 
dius in two weak and strong gravity region. It is seen 



Strong gravitational field in 



7 



that in the strong region (the close up part) the rele- 
vant analytical answer and the numerical solution are 
agree together while the analytical weak field approx- 
imation solution deviates from the numerical solution. 
The close up part of figures show that with increas- 
ing the radius and going to the weak filed region, the 
analytical solutions of strong filed approximation and 
numerical answers get separated from each other, and 
at last in the weak field region, i.e. x ^> f, the ana- 
lytical weak field answers coincide with the numerical 
solution, while the answers for the strong gravity region 
has a grate deviation from the numerical results in this 
region. 
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